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Skyrmions are topologically protected spin textures, characterized by a topological winding num¬ 
ber iV, that occur spontaneously in some magnetic materials. Recent experiments have demonstrated 
the capability to grow graphene on top Fe/Ir, a system that exhibits a two dimensional Skyrmion 
lattice. Here we show that a weak exchange coupling between the Dirac electrons in graphene 
and a two dimensional Skyrmion lattice with A/" = ±1 drives graphene into a quantum anomalous 

Hall phase, with a band-gap in bulk, a Chern number C — 2N and chiral edge states with perfect 

2 

quantization of conductance G = 2N ^. Our findings imply that the topological properties of the 
Skyrmion lattice can be imprinted in the Dirac electrons of graphene. 

PACS numbers: 


I. INTRODUCTION 


Graphene is a zero gap semiconductor in the brink of 
becoming a topological insulator. It is no coincidence 
that several predictions of topological phases in two di¬ 
mensional systems are based on some small modification 
of the model that actually describes graphene, namely, 
a tight-binding Hamiltonian for electrons in a honey¬ 
comb lattice at half filling. In a seminal papeilll, Hal¬ 
dane proposed that a honeycomb crystal with a suit¬ 
able magnetic flux decoration would become insulator 
with a quantized Hall response, without Landau levels. 
The ground-breaking proposal of the quantum Spin Hall 
phase by Kane and Mel^^ is also based on graphene 
with spin-orbit coupling, that is mathematically related 
to the Haldane model. A QSH-like phase was also pre¬ 
dicted for graphene under the influence of a perpendicu¬ 
lar magnetic field and ferromagnetic ordeJ^. In addition, 
the quantum Hall effect has been observed in graphene 
at a relatively low magnetic field at low temperature^, 
and even at room temperatur^ at high fields. Further¬ 
more, the combination of exchange interaction and spin- 
orbit coupling has been predicted to give rise to the QAH 
phas^. Here we propose a topological phase that, unlike 
all of the above, requires no spin-orbit coupling and no 
applied magnetic field. 

Topological insulating phases in two dimensions at¬ 
tract interest because of the very special transport prop¬ 
erties, such as the perfect conductance quantization of 
the quantum Hall phas^. These special transport prop¬ 
erties arise from the fact that topological phases have a 
gapped bulk and conducting edge states. Interestingly, 
QAH phases have chiral edge states for which intraedge 
backscattering is impossible. Thus, quantization of con¬ 
ductance of the QAH phase should be as perfect as the 
one observed in the quantum Hall effect. This contrasts 
with QSH edge states, for which only the total absence 
of time-reversal symmetry breaking impurities, such as 
nuclear spins or local moments, prevents backscattering. 
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FIG. 1: (a) Scheme of a graphene layer deposited on a 

skyrmion lattice (only the skyrmion cores are shown), (b) 
Skyrmions with A = ±1, that trigger a QAH. (c) Influence of 
the exchange to the skyrmions on the Dirac cones: spin split¬ 
ting and band-gap opening (d) Spin texture induced in the 
reciprocal space by the Skyrmion driven gap opening. Edge 
states between vacuum and 2d Skyrmion crystal (e), and be¬ 
tween two skyrmionic crystals with different chirality (f). 


As a result, the quanti zatio n of conductance in QSH and 
QSH-like phases is faJ^^^ from the perfection observed 
in graphene QH system^^. 

The notion that exchange interaction with the local 
moments with non-collinear spin textures affects severely 
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the kinetic energy of itinerant electrons, goes back 
to the proposal of the double-exchange mechanisnP^. 
Later, the role of no n trivial effects on the Berry phase 
was recognizecP^^, and the notion that non-collinear 
spin textures induce an effective orbital magnetic field 
that would lead to an anomalous Hall term was put 
forwarcP^H^. It was also shown by Ohgushi et 
that the double exchange model on a two dimensional 
Kagome lattice with non-collinear classical spins led to a 
quantum anomalous Hall phase with quantized Hall con¬ 
ductance. Recent experiments have demonstrated the 
possibility of growing graphene islands on top of a two 
dimensional skyrmion lattice hosted by a single atomic 
layer of Fe deposited on an Ir(lll) substrat^^. This 
experiment, together with the notion that non-collinear 
spin textures in duce very interesting effects on conduc¬ 
tion electrons IIH li 5 | 20t {22|^ motivate our study of graphene 
electrons coupled to a skyrmion lattice. Interestingly, we 
find that this system can exhibit the Quantum Anoma¬ 
lous Hall effect (Fig. le,lf), a topological phase of mat¬ 
ter that is being actively pursued in condensed matter 
physic^^ 

The paper is organized as follows. In Section H we 
present the tight binding model to study the skyrmion 
proximity effect in graphene. In section HI we show the 
non-trivial phases that arise by proximizing graphene to a 
triangular and rectangular skyrmion lattices respectively, 
in the weak coupling limit at half filling. Finally, in sec¬ 
tion IV we summarize our conclusions. 


Mi = Mofii are the electronic spin density operator and 
classical magnetization, respectively, evaluated at site i. 
The length of the magnetization field Mq is assumed to 
be the same for all sites. Unlike most of previous worlP^, 
we focus on the weak coupling limit J << t, adequate for 
proximity induced magnetism. Therefore, the exchange 
field act as a perturbation on the Dirac spectrum. 

The magnetization of a single skyrmion can be ex¬ 
pressed as a map from the plane described in polar co¬ 
ordinates r, (j) to the unit sphere, described in spherical 
coordinates 0, 


n = (sin 0 (r) cos ^( 0 ), sin 0 (r) sin ^( 0 ), cos 0 (r)) ( 3 ) 


where 0(r) is such that 0(r = 0) = 0 and 0(r > Rsky) = 
TT , where Rsky is the skyrmion radius. ^( 0 ) = N(j) + 7, 
where N is the skyrmion number that accounts for its 
vorticit}® and 7 is a phase that determines the helicity of 
the skyrmion. Since 7 can be gauged by a rotation along 
the z-axis, we take 7 = 0 without loss of generality. For 
the sake of simplicity, in the following we will assume a 
stepwise profile for the magnetization 


0hard(^) 


0 r = 0 
< 7 r /2 0 < r < Rsky 

TT r>Rskv 


( 4 ) 


We have verified that the results do not change quali¬ 
tatively if we use a smooth parametrization of the az¬ 
imuthal angle. 


II. MODEL 

Skyrmions are non-collinear spin textures (Fig. la and 
lb) characterized by a topological numbeil^ 



where n = ^ is the unit vector associated to the 

|M| 

skyrmion magnetization. Skyrmions are being very ac¬ 
tively studied in the context of spintronic^^^, and they 
have been found both in bulk compound^ and two di¬ 
mensional system^^. Here we propose a mechanism to 
induce the QAH phase in graphene, fully independent of 
the spin-orbit coupling of carbon. Namely, it is based on 
exchange interaction with a magnetically ordered surface 
that hosts a skyrmion lattice (Fig. la). 

Our starting point is the tight-binding model for elec¬ 
trons in graphene Rq = plus their exchange 

to an arbitrary magnetization field, that is treated clas¬ 
sically: 

n = no + Mi ( 2 ) 

i 

Here i labels the sites of the honeycomb lattice, J is 
the short-range exchange interaction constant. Si and 


III. QUANTUM ANOMALOUS PHASE AT 
WEAK COUPLING 

In this section we present the result for the weak ex¬ 
change limit (J << t). This situation is the one to be 
realistically obtained for graphene on top of a skyrmion 
lattice. We show that a topological gap opens for arbi¬ 
trarily small exchange coupling and independently on the 
type of skyrmion lattice, triangular or rectangular. 

A. Triangular skyrmion lattice 

We now consider the properties of graphene interact¬ 
ing with a triangular crystal of skyrmions. In Fig 2 b 
we show the energy bands for a 5 x 5 supercell with one 
skyrmion. It is apparent that exchange interaction with 
the skyrmions opens up a gap in graphene. This can be 
understood as follows. We write the exchange part of the 
Hamiltonian as 

V = +jY,Si-({M) + 5Mi') ( 5 ) 

I 

where (M) = ^ Mj is the average magnetization, 
and 5 Mi = Mi — (M) are the fluctuations. Ignoring the 
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FIG. 2: (a) Scheme of a triangular arrangement of skyrmions 
and (b) band structure of a 5x5 honeycomb supercell with a 
skyrmion with winding number equal to 1 and J = 0.3t. (c) 
Contour plot of the Berry curvature, with a small trigonal 
warping around the corners of the hexagonal Brillouin zone, 
(d) Positive Berry curvature localized at the region of band- 
gap opening, (e) Surface density of states on the termination 
of a semi-infinite graphene plane, showing two in-gap chiral 
states, (f) Anomalous Hall conductivity as a function of Fermi 
energy 


fluctuation term, the average magnetization induces a 
spin-splitting, shown in figure lb, where the bands have 
a well defined spin along the average magnetization. The 
resulting conduction band of one spin projection is de¬ 
generate with the valence band of the opposite spin pro¬ 
jection, defining a circle of degenerate points that, at half 
filling, happens to be the Fermi surface (see Fig. Ic). For 
non-collinear spin textures, 5M has terms orthogonal to 
(M), that open up the gap at the degeneracy circle, 

and creates a spin vorticity in the reciprocal space (see 
Fig. Id). A similar argument has been used by Qiao and 
coworkers in their proposaP for QAH in graphene with 
Rashba spin-orbit coupling and a exchange field. In our 
system the 5M^ part of the Hamiltonian plays the same 
role as the Rashba coupling in their model. However, 
we found that only skyrmions with topological charge 



FIG. 3: (a) Scheme of a rectangular skyrmion crystal, (b) 

Band structure of half filled graphene with a rectangular 
unit cell, coupled to the rectangular skyrmion lattice with 
J = 0.3t. (c) Berry curvature in the whole Brillouin zone, 
showing a same sign behavior which sums up to a C = +2 
total Ghern number, (d) Berry curvature along the k path 
shown in the band structure, showing a non-vanishing con¬ 
tribution in the anti-crossings points, (e) Surface DOS in a 
seminfinite geometry showing two chiral states, (f) Anoma¬ 
lous Hall conductivity as a function of Fermi energy 


N = are able to open a topological gap. 

The non-trivial nature of the gap of Fig. 2b can be 
anticipated from the Berry curvatur^^ profile shown in 
Fig 2d. A non vanishing Berry curvature arises at the 
crossing circle of the spin-split Dirac points (see Fig. 2c), 
where the spin-flip terms open-up a gap. 

The Hall conductivity c an be b e expressed as an inte¬ 
gral of the Berry curvatur^H^^ISlI; 

a^yiEp) = f{en{k),EF)^n{k)d?k (6) 


where f{en{k)) is the Fermi occupation function, and n 
labels the bands. When the Fermi energy, Ep, lies inside 
a gap, the Hall conductivity is proportional to a Chern 
number d^, which is an integer number C, resulting in 
a quantized Hall conductivity axy = Ce^ jh. Our calcu¬ 
lations show that when Ep lies in the gap opened by the 
coupling to the N = El skyrmions, the Chern number is 
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given by. 

C = 2N (7) 

This is the central result of this work: the topological 
winding number of the skyrmions is imprinted into the 
Dirac electrons. Several consequences follow. First, two 
chiral edge states are expected to occur at the boundaries 
of the crystal. This is confirmed by our calculations, us¬ 
ing a recursive Green function method^ to calculate the 
surface states, as shown in Fig. 2e. Second, an interface 
between two skyrmion lattices, with opposite skyrmion 
numbers N = 1 and N = —1, is expected to show 4 
interface states according to the index theorem (see fig. 
If). Third, when Ep lies inside the gap, the Hall con¬ 
ductance is quantizecPSl^ axy = C ^. Importantly, the 
Hall conductivity also takes large, but not-quantized val¬ 
ues, when effect Fermi energy lies close to, but outside, 
the gap (see Figs. 2f, 3f), on account of the finite Berry 
curvature integrated over the occupied state^^. 


B. Square skyrmion lattice 

We now address the question of the influence of 
the type of skyrmion lattice on the existence of the 
QAH phase. In particular, motivated by experimental 
result^, we consider a rectangular skyrmionic crystal, 
commensurate with the graphene lattice (see figure 3a). 
We find that the topological character of such system 
strongly depends on whether the valley mixing is an im¬ 
portant effect. 

The topological phase is observed for those unit cells 
in which the two valleys fold to different points in the 
Brillouin zone, which avoid intervalley mixing. As in the 
triangular case, the in-plane components open up a gap 
in the exchange split bands (Fig. 3b), which leads to 
a finite non-vanishing Berry curvature localized in the 
band crossing points (Fig. 3d). The Fermi surface in 
this case can be more complex than in the triangular 
skyrmion lattice due to the large folding of the reciprocal 
space, which is also reflected in Berry curvature Fig. 3c. 
The calculated Chern number is also C = 2N, as in the 
triangular lattice, which in a semi-infinite geometry gives 
rise to two copropagating edge branches (Fig. 3e). 

Away from half filling, when the Fermi energy no 
longer lies within the gap, the anomalous response is still 
non-vanishing up to energies 4 times the gap (Fig. 3f). 
The sign of the anomalous response is again the same for 
electrons and holes as in the triangular case. In compar¬ 
ison, upon entering the conduction (valence) band, the 
response can rapidly become negative and with a value 
close to 1, due to the presence of an edge branch located 
in the valence (conduction) band (see Fig. 3e), that co¬ 
exist with normal valence (conduction) states. 


C. Experimental verification 

The experimental verification of our proposal is not far 
from the state of the art. The optimal conditions to de¬ 
tect perfect edge transport associated to the QAH phase 
is induced by a skyrmion lattice hosted by an insulating 
magnetic material that couples to the graphene electrons. 
Three ingredients have been demonstrated in different 
systems. First, recent experimental results showing the 
possibility of growing graphene on the surface of Fe/Ir, 
a surface that hosts a skyrmion lattic^^. Second, and 
independent from the first, recent experiments show that 
non-quantized anomalous Hall is induced in graphene by 
proximity to a ferromagnetic insulating substrat^^. Fi¬ 
nally, skyrmion lattices have been observed in insulating 
chiral-lattice magnet Cu 20 Se 0 j^. 

The critical figure of merit to realize our proposal is 
the magnitude of the skyrmion induced gap. This is 
controlled by the strength of the exchange field of the 
underlying magnetic state. To gain some insight of the 
gap opening mechanism, it is convenient to change in¬ 
dependently the off-plane (exchange shift) and in-plane 
(spin mixing) components of the skyrmion texture. Our 
numerical calculations (see Fig. 4b) for the triangular 
lattice of hard core skyrmions {nz{r < Rsky) = 0) show 
that the gap satisfies A oc JinJojj where and 
are the magnitude of the in-plane and off-plane compo¬ 
nents of the exchange field (see Fig. 4b,4c,4d). In the 
case of constant exchange strength, a quadratic exchange 
dependence is obtained A oc J^/t. 

For reference, and given that there are no measure¬ 
ments of J, we take it from DFT calculations for graphene 
on top of BiFeO^ that give J = 70 meV, and t = 2.6 
eV. With the previous parameters, the topological gap 
would get a value of A 0.4 meV, within reach of trans¬ 
port spectroscopy. Taking these numbers, the window 
of Ep within which the non-quantized intrinsic Hall ef¬ 
fect would be sizable extends in a window of 1 — 3 meV 
around the Dirac energy. Our calculations indicate that 
the non-quantized Hall conductivity would be larger in 
the case of the square lattice (Fig. 3f). We also note that 
the sign of the skyrmion-induced Hall contribution would 
be the same for electrons and holes, in contrast with the 
conventional Hall effect. 


IV. CONCLUSION 

We have shown that, at half filling, graphene with a 
weak exchange coupling to a skyrmion lattice develops a 
quantum anomalous Hall phase, with gapped bulk and 
chiral edge states that should have perfect quantization. 
This occurs at least for two different skyrmion lattices, 
rectangular and triangular, and seems a generic feature as 
long as the skyrmion lattice does not produce valley mix¬ 
ing. The Chern invariant C that characterizes the QAH 
phase is given by the topological invariant that describes 
the individual skyrmions V, through the remarkable re- 
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J/t 




Joff/t Jin/t 


FIG. 4: (a) Scheme of a hard core skyrmion, showing an in¬ 

plane radial magnetic profile and an out-core off-plane mag¬ 
netism. Evolution of the topological gap at half hlling with 
the global exchange (b), only the out core exchange (c) and 
the core exchange (d). 


lation C = 2N^ valid for = ±1. Thus, graphene edges 
will have 2 chiral edge states, and graphene on top of an 
interface between two skyrmions lattices with opposite 
skyrmion number N = -\-l and N = —1 will have 4 chi¬ 
ral edge states. Our proposal is different from previous 
proposals of QAH phase in graphene because it requires 
no spin-orbit coupling and no magnetic field acting on 
the graphene electrons. 

Note added: During the final stages of the comple¬ 
tion of this manuscript, we became aware of a work pre¬ 
dicting a Quantum Anomalous phase for electrons on a 
square lattice coupled to skyrmion^^ in the strong cou¬ 
pling limit (J >> t). 
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FIG. 5: Graphene supercell over a triangular skyrmion crys¬ 
tal. Band structure (a,b) at hlling 1/4 in the large exchange 
limit (J=3t). As shown by the surface density of states (c) 
and the Berry curvature (d), the system develops a quantum 
anomalous Hall state with a total Ghern number C = 1 At 
large exchange, the half hlling state becomes trivial, whereas 
the hlling 1/4 becomes topological resembling the Haldane 
model. 


realistic in the case of graphene. In the strong coupling 
limit at half hlling, the system behaves as a topologically 
trivial magnetic insulator, where both valence and con¬ 
duction states of the majority spin are completely full. 
Since spin degeneracy is completely lifted in this case, the 
Fermi energy lies at the Dirac point at quarter hlling. In 
this case (Figs. [^, the low energy states are described by 
a gaped Dirac like spectrum and our calculations in this 
limit give C = 1, for A = 1. Therefore, this is topologi¬ 
cally different to the weak-coupling case discussed in the 
main text, and much closer to the original QAH phase 
proposed by Haldan^. 

The physical origin of this topologi cal ph ase can be 
understood as follows. The standar j^^FQH H] ro¬ 
tation is perf orme d, so that the exchange term is al¬ 
ways diagona^^^. The local nature of this transfor¬ 
mation gener ates a coupling to a gauge held in the hop¬ 
ping operatoj ^^FOtf H] that describes an effective inhomo¬ 
geneous magnetic held, that is responsible of the band- 
gap opening in this limit. In that limit, a skyrmion with 
topological number N creates an effective held equal to 
where = e/h is the magnetic hux quantunP. 


Appendix A: Strong coupling limit Appendix B: Trivial gap opening 


We discuss here the behavior of the system in the large 
J limit, where J » t and the spin-splitting of the bands 
is much larger than the bandwidth. This strong coupling 
limit has been considered in previous workJ^, but is not 


In this section we will discuss some of the situations in 
which a skyrmion texture will not give rise to a topolog¬ 
ical gap. First we emphasize that strong defects or inho¬ 
mogeneities both in graphene, the underlying skyrmion 
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FIG. 6: Unit cells (a,b) and band structures (c,d) for graphene 
over a fully in-plane triangular (a,c) and rectangular (b,d) 
skyrmion lattice. A gap is opened in the band structure, 
but the Chern number is identically zero, yielding a trivial 
insulator. 



lattice or the exchange coupling can give rise to a phase 
transition between the topological state to a trivial state. 
This problem is common to any engineered topological 
state, and its study relies on the particular features of 
the microscopic models, which depend on the underlying 
material that hosts the skyrmions. 

In the following we focus on homogeneous effects well 
captured by our phenomenological model that can ruin 
the perfect C = 2 state. We will discuss two situation, 
the case of a pure coplanar spin texture, and situations 
with strong intervalley mixing. 

1. Coplanar spin textures 

Here we consider the case of graphene coupled to copla¬ 
nar non-collinear spin textures, shown in Fig. We ad¬ 
dress the weak coupling limit (J << t) at half filling. In 
the case of a purely in-plane exchange field, the Hamilto¬ 
nian can be made real by a rotation onto the Pauli ma¬ 
trices (Jx^cFzi leading to a vanishing Hall response. We 
show in Fig. examples of band structures of triangu¬ 
lar and rectangular graphene unit cells subjected to fully 
in-plane exchange field. Although a band-gap opens, the 
calculated Berry curvature, and thereby the Chern num¬ 
ber vanish in both cases. 

The necessity of a non-coplanar spin texture can be 
easily understood in terms of the symmetry properties 
of the Chern number. Without loss of generality, the 
exchange field of a coplanar spin texture can be expressed 
in terms of the Pauli matrices cTx and cTz^ provided the 
exchange field is rotated to lie in the xz plane 

M ’ a = MxCFx + MzCFz (B1) 

turning the exchange term of the Hamiltonian purely 


FIG. 7: Unit cells (a,b), band structures (c,d) and Berry cur¬ 
vatures (e,f) for graphene over a triangular (a,c) and rectan¬ 
gular (b,d) skyrmion lattice. The commesuration of graphene 
with the skyrmion lattice make that both valleys are folded 
onto the F point. Even though a gap with non vanishing Berry 
curvature opens up, the sign changes along the Brillouin zone, 
and summing up all the contributions gives a vanishing Ghern 
number. In the present situation, the gap is dominated by in¬ 
tervalley mixing. 

real, and therefore also the full Hamiltonian. Since the 
Chern number is odd under conjugation 

K :C^ -C (B2) 

and a purely real Hamiltonian is invariant under conju¬ 
gation 

K = n (bs) 

the Chern number is identically zero (C = 0) for a co¬ 
planar spin texture. Thus, non-coplanar spin arrange¬ 
ments are necessary to induce the anomalous Hall phase. 

2. Intervalley mixing 

In some instances the skyrmion lattice will be com¬ 
mensurate with the carbon honeycomb lattice, the two 
valleys in band structure of graphene could be folded to 
a single point, allowing to intervalley mixing. In par¬ 
ticular, for the triangular lattice, the folding of the two 
valleys onto the F point takes place for 3n x 3n unit cells. 
For a rectangular unit cell, the folding takes place for a 
supercell 3n in the zigzag direction. In these situations. 
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intervalley scattering can open a trivial gap, as shown in 
Figl3 of such behavior. Interestingly, the Berry curva¬ 
ture can be non zero, but the Chern number vanishes. 


Importantly, this situation requires a fine-tuning of the 
skyrmion and carbon lattices. In general, this is not the 
case, and non-trivial gaps are to be expected. 
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